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ABSTRACT 

We investigate the Lya line transfer in nearby and high redshift starbursting galax- 
ies, where the effects of high optical depths and the role of dust in the scattering medium 
are expected to be conspicuous and should be treated in a very careful manner. We 
present our first results in a dustless, static, and uniform H I media with moderate Lya 
line center optical depths tq = 10'^"^. We assume the temperatures of the media to be 
T = 10^~^ K, and use a Monte Carlo technique. We investigate the basic processes of 
the line transfer and confirm the criterion of otq > 10^ for the validity of the diffusion 
approximation suggested by Neufeld in 1990, where o is the Voigt parameter. Adopt- 
ing the model for evolution of a galactic supershell suggested by Tenorio-Tagle et al., 
we performed calculations on the Lya line formation for each evolutionary stage of an 
expanding supershell. The emergent Lya profiles are characterized by the double peaks 
and the absorption trough at the line center. We found that the absorption troughs 
expected in most of the evolutionary stages are not wide enough to be observed with 
current instruments. However, the absorption trough in the Lya emission profile from 
an expanding recombining supershell can be marginally detected. 

Subject headings: line: formation — radiative transfer — galaxies: starburst — galaxy: 
formation 



1. Introduction 

With the advent of large ground-based and space telescopes, primeval galaxies at high redshifts 
(3 < z < 5) have been detected through a variety of methods. One of the most efficient detection 
methods is the Lyman break method which can be used to select high redshifted objects. Subsequent 
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spectroscopy revealed the existence of Lya features in either emission or absorption (Steidel et al. 
1996, 1999) in many of galaxies showing Lyman break. 

As well as being a redshift indicator, Lya bears information on the physical state of those 
galaxies. The emergent Lya profile is very sensitive to various factors including the kinematics, 
scattering optical depth, and dust content. Hence, it is important to investigate the Lya line 
formation for its cosmological applications. 

The Lya profiles of many Lyman break galaxies in the early Universe may be classified into 
three types : (1) symmetric profile, (2) asymmetric or P Cyg type profile, (3) broad absorption 
profile extending to the Lorentzian wings. It is known that the Lya profiles of nearby starbursting 
galaxies also fall into one of the above three types (Kunth et al. 1998). It seems to be a general con- 
sensus that the P Cyg type Lya profile is attributed to the expanding neutral medium surrounding 
a Lya source (Lee & Ahn 1998; Kunth et al. 1998; Ahn 2000). This is supported by the existence of 
the interstellar absorption from low ionization species such as O I A1302 and C II A1335 which are 
blueshifted relative to the Lya line center (Kunth et al. 1998; Spinrad et al. 1999). It is also proba- 
ble that the profile is formed in galactic superwinds. Assuming that a galactic superwind possesses 
a Hubble type outflow, Lee k, Ahn (1998) calculated the emergent Lya profiles and investigated 
their spectropolarimetric properties. 

The P Cyg type Lya profile can also originate from a supershell which is a remnant of strong ex- 
plosions in a central star-forming cluster. Prom these pieces of observational evidence, Tenorio-Tagle 
et al. (1999) proposed a model for the evolution of galactic supershells based upon hydrodynamical 
simulations and ionization calculations. However, they presented only a schematic description of 
the Lya profiles in each evolutionary stage, so more systematic and quantitative investigations are 
desired. 

Lya is one of the strongest resonance lines found in various emission line objects and therefore 
has been the subject of intensive research for a long time. Both numerical and analytical studies 
have been carried out by a number of authors. Recent advances in computer technology have made 
it possible to explore the parameter space that was not previously accessible. 

The main objective of this study is therefore to understand the Lya line formation in a thick 
medium by using a sophisticated Monte Carlo method, and to apply the results to the models for 
the evolution of superbubbles in starbursting galaxies. This is the first paper in a series in which 
the Lya line transfer in an optically thick medium is investigated. In this paper, we consider only 
moderately thick and dustless media, and defer the cases of extremely thick media with and without 
dust to subsequent papers. 

This paper is composed as follows. In Section 2, we describe the model for starbursting galaxies 
and our Monte Carlo code. In Section 3, we show our results, and in the following section we discuss 
the astronomical applications of our works. The final section summarizes our major findings. 
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2. Model and Method 
2.1. Model 

Throughout this paper we consider a uniform, static, and plane parallel medium consisting of 
pure neutral hydrogen atoms. We characterize the thickness of the medium by the Lya line center 
optical depth denoted by tq from the midplane to either side in the z-direction. We assume that a 
young and massive star cluster is located at the center of the slab. This model describes the young 
starbursting galaxies detected by the Lyman break method. 

Recently Tenorio-Taglc et al. (1999) modeled the starbursting galaxies by a young star-forming 
cluster surrounded by a supershell, which evolves with age. According to their Fig. 8, we can 
summarize the evolution of superbubbles in starbursting galaxies as follows. 

Stage (a) - In the initial stage of starbursting galaxies, the surrounding medium is neutral and static 
with respect to the central photon source. During this phase {tage < 1-5 Myrs), the ultraviolet 
photons from the central cluster ionize the surrounding neutral medium. The supernova explosions 
also deposit the kinetic energy to the medium surrounding the star cluster. Subsequently an 
expanding shell, which captures the UV photons inside the shell, develops. 

Stage (b) - The shell breaks up by the Rayleigh- Taylor instability, and the UV photons can freely 
ionize the surrounding medium. Its ionization front usually gets to the outer edge in tage — 
2.5 Myrs. At this time, a conical H II region of partially ionized gas is formed. The ionized 
medium starts to expand due to the increase of the internal pressure and the mechanical energy 
injected by the central star cluster, thereby forming a supershell. 

Stage (c) - The supershell becomes radiative and emits Lyct photons. The neutral column density 

of the expanding shell increases and its expansion speed decreases. 

Stage (d) - Finally the neutral supershell forms, and the expansion almost ceases. 

In view of the above scenario, the starbursting galaxies are composed of the central Lya 
emitter or a young star-forming cluster and the surrounding media. There are three ingredients 
in modeling starbursting galaxies. The first is the hydrodynamical evolution of a supershell, the 
second is the photoionization, and the last is the state of the photon source. The former two 
ingredients are reflected in Tenorio-Tagle et al. (1999), and the last one opens two possibilities in 
the input spectrum. When the major contributers of UV photons (i.e., O, B stars) arc still alive, 
the input spectrum consists of the Lya emission and the UV continuum. Therefore, in order to 
model the Lya line formation in each evolutionary stage of the starbursting galaxies, we consider 
following six types of configurations. 

Case I - the scattering medium is neutral and cold (T = 10— 100 K), and the Lya source is located 
at the center of the slab. 

Case II - the scattering medium is partially ionized and hot (T = 10^ K), and the Lya source is 
located at the center of the slab. 

Case III - the scattering medium is partially ionized and also acts as a Lya emission source. 
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Case IV - a continuum source is located at the center of the scattering medium which is partiahy 
ionized and has a temperature T = 10^ K. 

Case V - a continuum source is located at the center of the scattering medium which is cold with 
temperature T = 10 K. 

Case VI - there is a Lya source at the center and the scattering medium is an expanding neutral 
shell with T = 10 K. 

Stage (a) is described by Case I and Case V, because the source of Lja and UV continuum 
is surrounded by the static, cold (T = 10 K), and extremely thick (tq = 10^~^) neutral medium. 
Unfortunately, our code becomes very slow when aro > 10^, where a is the Voigt parameter. We 
need a special scheme to overcome this limitation. In the present study, we consider cases with 
To = 10^-6 only. 

Assuming that the Lya photons in Stage (b) arc generated in the central star-forming region 
that is partially ionized, this can be modeled by Case II and Case IV. 

Stage (c) has two Lya emission sources. The major one is the central star cluster, and the 
minor one is the expanding radiative supershell. The major emission source is similar to that in 
Stage (b), and thus can be modeled by Case II and Case IV before the H I column density of 
the supershell becomes large enough to produce the P Cyg type absorption blueward of the major 
peak. The conical H II region of Stage (b) can be described by Case III. Case III also can describe 
the line formation in H II regions where young stars and interstellar media are all mixed. 

In the earlier epoch of Stage (d), the expansion speed of the neutral supershell is large, and 
there appears a P Cyg type Lya emission line. This era can be described by Case VI. Case V 
models the latter epoch of Stage (d) when the Lya sources or O stars have been extinguished. 
However, extremely thick cases will be briefly treated in the current work and more details will be 
presented in the future (see also Lee Sz Ahn (1998) for the galactic superwind). 



2.2. Previous Studies of Lya Transfer 

Since we arc interested in the emergent Lya profiles in various situations, we first introduce 
the dimensionless parameter x defined by 

X = Au/AuD = {y- yo)/AuD, (1) 

which describes the frequency shift from the line center in units of the Doppler shift Auu = 
i^o{vth/c). Here vth is the thermal speed of the scattering medium and c is the speed of light. 

The Lya optical depth for a given system is 

Tx = ToH{x,a), (2) 

where H{a,x) is the Voigt function or Hjerting function, a = 4.71 x 10~^T^^^'^ is the Voigt param- 
eter, and T4 represents the temperature of medium in units of 10^ K. The line center optical depth 
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is related to the H I column density Nhi via 



TO = 1.41 Tl 



1/2 



Nhi 



(3) 



10i3cm-2 



Since the normalized profile function (l){x) is related to the Voigt function by 



^(x) = —=H(x,d), 



(4) 



the line center optical depth used by previous authors (tq) is related with that defined in this paper 
(to) by To = ^To. 

Theoretical studies on the Lya line transfer in an optically thick and static medium have a long 
history. Unno (1955) formulated the Lya line transfer problem, and Osterbrock (1962) proposed 
a simple physical picture for understanding the resonance Une transfer in a thick medium. Adams 
(1972) revised Osterbrock's picture and gave a detailed explanation of the transfer mechanism. 
Noting that each scattering event is accompanied by an average shift in frequency of order the 
thermal Doppler width Az/£), he emphasized the diffusive nature of the transfer process both in 
frequency space and in real space introducing the terms 'single longest flight' and 'single longest 
excursion.' 

Neufcld (1990) divided the optical depths into three regimes: the slightly thick optical depth 
< 10'^, the moderately thick optical depth 10"^ < To < 10^/a, and the extremely thick optical 
depth To > 10^ /o. Here the range of a is between 10~^ and 10~^ for a scattering medium with 
10 < T < 10^ K. 

Adopting the diffusion approximation, an analytic solution for the extremely thick case was 
given by Harrington (1973), and a more general solution was obtained by Neufeld (1990). In 
moderately thick media, however, it is difficult to solve the problem analytically, because the 
diffusion approximation is not valid in this regime. Hence, a Monte Carlo method is useful to solve 
the problem of the Lya transfer in moderately thick media. 

Avery & House (1968) investigated the resonance-line scattering by a Monte Carlo method. 
In their study, they focused on the difference between the complete redistribution and the partial 
redistribution. Their calculations reach the condition To = 10^ and a = 4.3 x 10^'^, which turned out 
to be about the critical optical depth between the extremely thick and the moderately thick cases. 
Gould & Weinberg (1996) also investigated a similar radiative transfer problem. However, they 
concentrated on the observational feasibility including extinction in a partially ionized envelope 
of Lyman limit galaxies, and the detailed line formation mechanism was not fully described. In 
addition to these works, there have been several investigations using a Monte Carlo method for the 
resonance line transfer in an optically thick and static medium in the literature (e.g. Meier &; Lee 



1981). 
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2.3. The Monte Carlo Technique 

In this study, we investigate the problem by means of a Monte Carlo technique. A detailed 
description of this code was presented in the previous paper (Ahn et al. 2000). Compared with 
those of other previous investigators, the major advantage of our code is that it incorporates all the 
quantum mechanics associated with both resonant and wing scatterings of Lya. Here, we neglect 
any collisional broadening in the upper state, and assume that the frequency of the scattered photon 
is identical to that of the incident photon in the rest frame of the scattering atom. This is Case 
II frequency redistribution in the atom's frame introduced by Hummer (1962) (see also Mihalas 
1978). Hence, in the code we transform the necessary quantities between the observer's rest frame 
and the scattering atom's rest frame. 

According to Stenflo (1980), the scattering phase function associated with wing scattering 
of Lya is identical with the classical Rayleigh scattering which yields 100 % polarization for 90° 

scattering. However, the resonant scattering of Lya associated with the transition lsi/2 — 
is characterized by the isotropic scattering phase function yielding totally unpolarized scattered 
radiation, whereas for the lsi/2 ~ 2^3/2 transition the scattering phase function is characterized 
by the maximum polarization Pmax = 3/7 for 90° scattering of unpolarized incident photons. 
Therefore, in our Monte Carlo code, we computed the probabilities for these scattering types at 
each scattering event. 

The absorption profile can be written as the oscillator strength-weighted sum of two Voigt 
functions, 

where 

and vi and 1^2 are the resonance frequencies corresponding to the lsi/2 — and lsi/2 ~ '^P3/2 

transitions. 

The scattering type can be determined from the scatterer's local velocity, of which the distri- 
bution f{y) is given by up to normalization 

f^y^ = ^^1^ (a^ + ( J - yr + + («2 - y?) " 

Here y = v/vfh is the velocity in units of the thermal velocity, and A^^ is the normalization constant. 
From the definition of the Voigt function, we see that 



= 3[H{a, Ul) + 2H{a, U2)]~\ 



(8) 
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We may loosely say that the scattering is resonant with the transition lsi/2 when 
in the rest frame of the scatterer the incident photon has the wavelength within several natural 
widths, say 10a; i.e., \ui — y\ < 10a. Therefore, the probability for a lsi/2 ~^ '^Pi/2 resonance 
transition is given by 

rui+lOa 

Pri - / dyfiy) 

J ui — lOa 

~ e-^'/[H{a,ui) + 2H{a,U2)], (9) 

because the change in the exponential part is negligible in the integration interval. A similar 
consideration gives the probability for lsi/2 '^P3/2 transition, 

rU2+Wa 

Pr2 ^ / dyf{y) 

J «2 — lOa 

~ 2e-"2/[i/(a,ni) + 2i/(a,U2)]. (10) 
Therefore, the probability that a given photon is scattered in the damping wings is 

Pnr = 1 — Pri — Pr2- (11) 

In the code we determine the scattering type in accordance with the probabilities Pri,Pr2 and 
Pnr- If scattering is chosen to be resonant, then we set u = x. Otherwise, the scattering occurs in 
the damping wings, and u is chosen in accordance with the velocity probability distribution given 
by Eq. (7). 

We assign the propagation direction kf of a scattered photon in accordance with one of the three 
phase functions. The scattered velocity component vj_ perpendicular to the initial direction ki on 
the plane spanned by fcj and kf is also governed by the Maxwell-Boltzmann velocity distribution, 
which is numerically obtained using the subroutine gasdevQ suggested by Press et al. (1989). 
The contribution Ax of the perpendicular velocity component v± to the frequency shift along the 
direction of kf is obviously 

Ax='-^[l-{ki-kff]'/\ (12) 

Therefore, the frequency shift X f of the scattered photon is given by 

Xf = Xi - u + u{ki ■ kf) + ^[1 - [ki ■ fe/)Y^^ (13) 

where Xi is the frequency shift of the incident photon. 

In each scattering event, the position of scattered photon is checked. If the photon escapes from 
the medium or \z\ < tq, we collect that photon according to its frequency and escaping direction. 
The frequency bin is chosen to be 0.25 < Ax < 2 considering the computing speed and the spectral 
resolution required to describe the details of the emergent profiles. The whole procedure is repeated 
to collect typically about 10^ photons in each frequency bin. This code, being very faithful to the 
quantum mechanics of the scattering of Lya, we may compute the polarization of Lya. This will 
be presented in the future (see also Lee & Ahn 1998). 
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3. Line Formation Mechanism 

3.1. Emergent Profiles 

In this section, in order to test our Monte Carlo code, we compute emergent profiles for the 
thick cases with various optical depths, and compare them with the analytic solutions considered 
by Harrington (1973) and Ncufcld (1990). They investigated the diffusive process of the resonance 
line transfer in an plane-parallel, extremely thick, and static medium with a Lya source in the 
central plane. We performed Monte Carlo calculations for the same cases, and compared results 
with the analytic solutions. Prom the comparison, we can determine the critical optical depth that 
divides the moderately and the extremely thick optical depths, where in the latter case the diffusion 
approximation is can be applied. 

By adopting the Eddington approximation, Harrington (1973) introduced the diffusion equa- 
tion for the angle-averaged intensity J(f, x) 

^ + ^ = -3^(.)^, (14) 

where E{f, x) is the photon generation rate (per unit mean optical depth per unit Doppler width 
per unit area) and the frequency parameter a is defined by 

(15) 



a = / dx'/(t){x'). 

Jo 



For the case of a monochromatic source at the line center located at the midplane of the slab 
with E{f,x) = y/3/2S{f)S{a) and (afo)^/'^ > 1, Harrington (1973); Neufeld (1990) derived the 
following analytic solution for the emergent mean intensity J, 

^ = \ (^(^\ 

^ ^"'""^ 24 afocosh[(7rV54)V2(|a;3|/a^Q)]- ^ 

We performed Monte Carlo calculations for the same case, and Fig. 1 shows our results. Notice 
that we omit the profiles in the blue part because the profiles arc symmetric with respect to the 
origin (x = 0). The solid lines represent the results of our Monte Carlo calculations, and the 
dotted lines represent the analytic results obtained by Neufeld (1990). Our results for (tq = 10^, 
a = 4.71 X 10~^) and (tq = 10^, a = 1.49 x 10~^) are in good agreement with the analytic solutions. 
However, the result for the case (tq = 10^, a = 4.71 x 10~^) shows noticeable discrepancy. Our 
Monte Carlo solution appears to be translated from the analytic sohition toward larger x by an 
amount of Ax ~ 1, where wing scatterings become important and (p{x) can be well approximated 

by 

<Pix) c ^. (17) 
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In the cases of dust free media, no loss of Lya line photons is permitted so that the flux is 
conserved. The profile function (j){x) can not be approximated by Eq. (17) at x < 3, but more 
adequately described by ~ (a/7r)/(l + x^). Therefore, Neufeld's calculation underestimates 
the number of core photons that are removed and ultimately redistributed to the wing regimes. 
This reduces the number of diffusively transferred wing photons, and hence the emergent spectra 
show significant discrepancy for moderately thick media. 

According to Neufeld (1990), the diffusion approximation, Eq. (17), is valid only for afo > 10^. 
For the former two cases, we have afo = 562, 264 for (tq = 10^, a = 4.71 x 10~^) and (tq = 
10^, a = 1.49 X 10-2), respectively. On the other hand, afo = 56.2 for (tq = 10^, a = 4.71 x 10"^). 
Therefore Neufeld's condition for the extremity of optical depth is marginally valid in the former 
two cases, and this also confirms that our code operates very well. 

3.2. Last Scattering Positions 

The last scattering positions of escaping photons are investigated in order to understand the 
line transfer processes in detail. We denote by the z-coordinate values of the last scattering 
positions of escaping photons and plot Pz/tq in Fig. 2. We identify two kinds of peaks in the 
distribution, one at the central region (Pz/tq ~ 0), and the other located near the boundaries 
{\Pz\/to ~ 1). We also find that the number of photons escaping from the center decreases as tq 
gets larger and the scattering medium colder. 

According to Adams (1972), in a slightly thick medium (tq < 10^), an initial core photon 
(xo = 0) experiences a large number of core scatterings until being scattered at the wing frequency 
by an atom that moves much faster than the thermal speed of the medium. When a wing scattering 
happens, the photon is ready to escape from the medium, and the escape is achieved by 'a single 
longest flight'. During the core scatterings, the mean free path of the photon is so small that the 
photon can not travel much. Hence, these photons constitute the central peak in Fig. 2. 

On the other hand, in a moderately thick medium (10^ < tq < 10^/a), a wing photon can 
not escape from the medium by a single wing scattering, but it can escape after only a few wing 
scatterings. However, during these few successive wing scatterings, the photon may become a core 
photon again by the so-called 'restoring force,' and then it experiences a large number of core 
scatterings. While these processes are repeated, a photon wanders about the medium before it 
escapes. By a Monte Carlo technique, Ahn et al. (2000) showed that this process, which they 
named 'a wandering,' indeed operates. In this regime, the last scattering position of escaping 
photons is distributed rather uniformly, as shown in Fig. 2. 

In the case of extremely thick media, Lya photons escape by the transfer process similar to 
that for the moderately thick medium case. They alternatively experience a large number of core 
scatterings and a series of wing scatterings. However, when the medium is extremely thick and the 
temperature of the scatterers is low, wing scatterings become more important in the line transfer 
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and the number of successive wing scatterings becomes large. Hence, the escape of photons is 
mostly achieved during the series of wing scatterings. Adams (1972) called this process 'excursions' 
and their escape is achieved by 'a single longest excursion'. 

As a result, photons gradually smear out near the boundary of the medium, and they escape 
by a single longest excursion. These photons form the peak at |Pz|/to — 0.95 and the small portion 
in 0.2 < |P2;|/to < 0.8 in Fig. 2. In view of these excursion processes, we can understand that the 
central peak is suppressed and the outer peaks dominate, as tq of the medium gets larger and the 
temperature of the scattering medium gets lower. 

3.3. Last Scattering Path Length 

Next, we investigate the distribution of photon path lengths just before escape. In Fig. 3 we 
show our results, where the horizontal axis represents the path length r/ which is the optical depth 
traversed by an escaping photon in the slab after its last scattering in units of tq. In Fig. 3 one 
may notice that there are three kinds of peaks. The first peak appears at r/ « 0, the second at 
1 < Tf < 2, and the third at r/ ~ 2. 

We may argue that the first peak is formed by the photons that have experienced 'single 
longest excursions'. Considering the processes mentioned in the previous subsection, these photons 
have traveled in real space through excursions and therefore escape happens preferentially for those 
photons that have diffused spatially near the boundary of the scattering medium. 

The second peak is due to the photons that escape by the 'single longest flight' from the central 
region of the medium. We may explain the shape of the distribution curve for this peak as follows. 
Since a large number of core scatterings isotropize the radiation field, as we will show in Fig. 4, the 
angular distribution of the emergent photons is approximately given by -P(/it) oc /x in the vicinity 
of the photon source. Here, ^ is the cosine of the angle between the outgoing wavevector and the 
normal direction of the scattering medium. Therefore, the escape by a single longest flight can 
happen in the direction /j, when the path length satisfies r/ ~ to/|//|. This implies that the path 
length distribution /(r/ = to/|/x|) oc \fi\ oc tq/t/ for r/ > tq, where /x G [0, 1]. In addition, a minor 
contribution to the second peak is due to the photons last scattered in the intermediate region, 
0.2 < IP^I/ro < 0.8. 

The third peak is contributed by the back-scattered photons. The frequency diffusion process 
drives these back-scattered photons to acquire very large wing frequencies, at which the medium 
becomes transparent enough for them to escape readily. Hence, this peak appears only when the 
wing scatterings become important. Therefore, as tq gets larger and the scattering medium colder, 
the third peak gets stronger and the first and second peaks become weaker. In particular, when 
To = 10^, T = 10 K, no second peak appears. 
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3.4. Dependence of the Mean Scattering Number on the Initial Frequencies 

In this subsection, we investigate the dependence of the mean number of scatterings {N) on 
the initial frequency xq and the Voigt parameter a. This was investigated by Avery & House 
(1968) using a weighting scheme. They introduced the weighting scheme in order not to count 
those photons that arc initially in the wing regime and directly escape without being scattered, 
and hence to prevent the {N) distribution from varying abruptly at 0. However, we did not adopt 
this scheme, and show our results in Fig. 5. Here we note that the areas under the curves are 
normalized to the same value because no dust is present in the medium and the weighting scheme 
is not used. 

Although, at small {N), the emergent flux appears rather fluctuating, our results are in quali- 
tative agreement with those obtained by Avery & House (1968). We see that the number of directly 
escaping photons increases as the initial photon frequency gets larger and the medium gets colder. 

We explain these facts as follows. The initial wing photons with xq > 2.5 experience coherent 
wing scatterings and can escape more easily with a small (N). However, the so-called 'restoring 
force' renders some of wing photons to be core photons, which will subsequently experience a large 
number of core scatterings. Therefore, we have two local maxima in (A^) of escaping photons and 
we get less photons with large {N) when xq is large. The larger fraction of directly escaping photons 
in the cold medium means that the smaller width of the Doppler core and the larger probability of 
wing scatterings in the Voigt profile play a main role in the line transfer. 

4. Observational Ramifications 

In this section we apply our code to model the Lycc line formation in starbursting galaxies. 
We adopt the models described in Section 2, and show our results of the Monte Carlo calculations. 

4.1. Case I - Hot Source and Cold Scatterers 

In Case 1, a hot source is embedded in a cold scattering medium with the line center optical 
depth To = 10^~^. Here, two types of sources are considered. The first is the O, B stars in a young 

star cluster, and the other is the central region ionized by them. However, for both sources the input 
profile can be assumed to have a width of about 10 km s^^. This is because a typical young cluster 
has a one dimensional velocity dispersion of 10 < ct^ < 20 km s~^, and the typical temperature 
of the ionized region is T = 10^ K or cr^ 10 km s~^. Due to the lack of our knowledge about 
the source, we just assume a Gaussian profile (f){x) oc exp(— x^/x^), for the Lya emission from the 
source, where Xg ~ 30. The value of Xs was chosen such that the line width of source is greater 
than the Doppler width of the scattering medium (T = 10 K) by a factor of 10^/^. We have not 
considered the microturbulence. 
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We show our Monte Carlo results in Fig. 6. The emergent profiles are characterized by double 
peaks and a central absorption trough, which are superposed on the broad source profile. The full 
width of the absorption trough is found to be Ax ~ 6, which corresponds to AA ~ 0.01 A in the 
rest frame. The value of AA is weakly dependent upon the temperature of the scattering medium. 

For those objects located at z = 3, the width of the absorption trough becomes 0.04 A. 
Considering that the spectral resolution of the Low Resolution Imaging Spectroscope (LRIS) on 
Keck I is about 2— 5 A (Oke 1995), the detection of the absorption trough profile may be impossible 
with LRIS, but it will be possible in the near future. 

The total hydrogen column density of normal spiral galaxies and dwarf galaxies ranges from 
Nhi = 10^9~22 (>j^-2^ which corresponds to the Lya line center optical depth of tq = 10^~^. Hence 
for such realistic broader absorption trough is expected, which can be more readily observ- 

able. Moreover dust can play an important role in widening the absorption trough. However, the 
Monte Carlo method becomes highly inefficient for the extremely thick cases with aro > 10'^, be- 
cause the number of core scatterings per excursion increase enormously. More detailed investigation 
on this issue will be reported in forthcoming papers. 

The two peaks around the central frequency are a consequence of the assumption of the dustless 
medium. The photons constituting the peaks have traversed very long distances. In this case, they 
are prone to be destroyed by even a small amount of dust that may co-exist in the scattering medium 
(Chen & Neufeld 1994). Hence we may make use of their existence and strength to measure the 
degree of metal contamination or the dust contents of the medium. We will report the results on 
this topic in the future. 

4.2. Case II - Monochromatic Source and Hot Scatterers 

When the scattering medium is partially ionized and the central Lya source is the O, B stars 
or the central dense H II region, we may safely assume that the source is monochromatic. This is 
because the width of the source is similar to the Doppler width of the scattering medium. Hence, we 
put a monochromatic and isotropic source with xq = at the center of the plane-parallel medium, 
where the temperature of the scattering medium is set to be T = 10^ K, typical value for H II 
regions. 

Our results are depicted in Fig. 7, where we show the effects of various line center optical 

depths. As tq gets larger, the peaks move farther away from the line center. When tq < 10^, the 
black absorption at the line center disappears due to the slab geometry whose covering factor is 
unity. We note that the moderate optical depth regime starts from tq > 10^, which can be seen 
from the appearance of the black troughs in Fig. 7. 

We show the effect of the Voigt parameter on the line formation in Fig. 8, where the case for 
To = 10^ is shown in the lower panel, and the case for tq = 10^ in the upper panel. Here, since 
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a = 4.71 X 10-'^(T/10'^K)-y^, a = 4.71 x 10"^ corresponds to T = 10^ K, a = 4.71 x 10"=^ to 
T = 10^ K, and a = 1.49 x 10~^ to T = 10 K, which are representative of the three phases of 
interstellar medium. 

In the figure, we see that the peak gets broader and the peak center moves farther from the line 
center as the scattering medium gets colder. This is because the radiative transfer is more affected 
by wing scatterings as a gets larger. We note that the increase of a affects the line formation in 
a similar way as the increase in tq in the extremely thick cases, where wing scatterings dominate 
the line transfer. We also note that the broad line spectra can be seen only for examples where the 
parameters approximately satisfy the condition (aro)^/^ S> 1, which is described in Section 3.1 and 
also discussed by Neufeld (1990). 

In this case. Ax = 1 corresponds to AA fa 0.05 A in the emitter rest frame. Therefore, the 
absorption trough has a width of AA ~ 0.3 A, which is expected in the ultraviolet spectra of nearby 
starbursting galaxies. When those objects arc located at z = 3, the width of the absorption feature 
is about 1.2 A which can be marginally detected by the Keck LRIS. 

Tenorio-Tagle et al. (1999) predicted only the symmetric single peak in Stage (b). We found 
from our calculations that Lya may have double peaks if the density of the central H II region is 
high enough to form a partially ionized medium. However, the validity of this prediction is largely 
dependent upon the existence of partially ionized zone in the star cluster. This may be checked by 
exploiting hydrodynamical simulations and calculating the ionization structure. We also note that 
dust may seriously alter these features. However, these topics are beyond the scope of this paper, 
and we concentrate on the radiative transfer problem in a pure hydrogen medium. 

4.3. Case III - Uniformly Distributed Source 

We also performed calculations for the case of uniformly distributed sources in the slab. This 
case can be applied to the recombining shell at the epoch (c) in Fig. 8 of Tenorio-Tagle et al. (1999). 
We distribute Lya sources at {Px,Py,Pz) = (0,0, r^), where Tg = {2R — l)ro and i? is a uniform 
random number in the interval [0, 1]. The number of sources are determined by the following two 
constraints. The uniformity of the sources should be well described by the Monte Carlo code, and 
the S/N ratio of the emergent profiles should be high enough to give us meaningful results. In our 
Monte Carlo computations approximately a few hundred photons per frequency bin with Aa; = 0.25 
are found to satisfy these constraints. 

Each source radiates monochromatic Lya photons isotropically. Since the scattering medium 
is partially ionized, the temperature of the medium is set to be T = 10"^ K. Moreover, since the 
Lyman limit opacity is smaller than the Lya line center optical depth by a factor of 10^""^, we 
performed the Monte Carlo calculations for this case by setting the Lya line center optical depth 
To < 10^. In other words, we assume that this case can be applied only before the recombination 
shell becomes extremely thick. 
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In Fig. 9 we show Lya profiles for 10 < tq < 10^ and temperature T = 10^ K. In the figure we 
can see the effects of the hne center optical depth on the emergent profiles. When compared with 
the results for the cases of the midplane source, no apparent discrepancy was found for tq = 10^. 
For 10 < To < 10^, the inner part of the profile gets broader, and the number of photons near the 
line center increases. This is caused by the contribution of photons at the shallow region, which can 
escape more easily without much frequency diffusion. However, when the line center optical depth 
is large, frequency diffusion takes place, and the emergent profiles become similar to the midplane 
cases. 

In Fig. 10, we also show the effect of the Voigt parameter on the line formation for the 
distributed source. We can see that the emergent profiles for the distributed source have a stiffer 
inner part of the peak than those for the midplane source. This trend is more conspicuous in the 
case of To = 10"^ than tq = 10^. In the case of tq = 10^, the escape of photons takes place mostly by 
single longest flights, and is less affected by the temperature of the scattering medium. However, 
when To = 10^ and T < 100 K, wing scatterings become dominant in the line transfer. Hence, 
photons smear out near the surface of the medium, and escape easily without much frequency 
diffusion. Therefore, the peak profile shows a stiffer inner part and the maximum flux occurs closer 
to the line center than in the midplane source case. 



4.4. Case IV &; V - Continuum Source 

When there are O, B stars alive in the central star- forming region, their continua also contribute 
to the Lya source photons, in addition to the Lya recombination line photons. Thus far only the 
line source has been considered. Now we consider the UV continuum as a photon source. Then 
the emergent profile is the sum of these two components, whose ratio is dependent upon the stellar 
evolution history of the central cluster. 

In Case IV, we assume that a continuum source with a flat spectrum is located at the center 
of the slab, and that the scattering medium is partially ionized or hot (T = 10^ K). This case is 
meaningful only when the total line center optical depth tq < 10^, and therefore in this study we 
consider 10 < tq < 10^ in order to safely cover the meaningful range of tq. 

In Case V, we assume a flat UV continuum as the photon source and a cold (T = 10 K) 
medium as the scatterer. In this case, the scattering medium is expected to have the line center 
optical depth tq = 10^~^, but due to the limitation in computing speed of our present code, we 
perform calculations only for tq = 10^~^. 

The results for Case IV and V are depicted in Fig. 11 and Fig. 12, respectively. For both cases 
IV and V, there appear double peaks with a broad absorption trough at the line center. For Case 
IV, the width of the trough is about 0.3 — 0.5 A. When the objects are located at the redshift z = 3, 
the width of the absorption trough is about 1.2 — 2.0 A, which is marginally observable with the 
large ground-based telescopes. The spectral features obtained in Case V are qualitatively similar 
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to those in Case IV, except for the smaller widths 0.01 — 0.03 A. When the source is located at 

the redshift z = 3, the width of the absorption trough ranges from 0.04 — 0.13 A. We expect that 
the absorption trough becomes mucli wider when we consider the effect of dust absorption and the 
high H I column density which is more realistic than the models considered in this work. 

5. Summary and Discussion 

In order to predict the Lya profiles from starbursting galaxies, we investigated the Lya line 
transfer in a moderately thick and dustless medium using a Monte Carlo code. Our code is faithful 
to the quantum mechanics associated with the resonant and wing scattering of Lya, and the partial 
frequency redistribution is accurately treated. 

We compared our results with the analytic solutions obtained by previous investigators, and 
found that the results are in excellent agreement. We confirmed the line transfer mechanism such 
as a single longest flight, a wandering, and a single longest excursion as proposed by Adams (1972). 
Comparing our emergent profiles with the analytical solution of Neufeld (1990), we found that the 
wing approximation does not hold for moderately thick media 10^ < tq < 10^/a, and that it is 
valid only for the regime otq > 10^. 

Next, we applied our results to the Lya line formation in starbursting galaxies, in relation with 
the model given by Tenorio-Tagle et al. (1999). The following is the summary of the emergent Lya 
profiles for each evolutionary stage of a superbubble in starbursting galaxies. A brief summary is 
also presented in Table 1, and the spectral capabilities of typical instruments are shown in Table 2 
for comparison. 

(a) In the early phase of a starbursting galaxy, the galaxy has a static and cold (T = 10 K) 
envelope, and the line width of the source can be set to be 10 km s^^. We obtained the emergent 
profiles consisting of double emission peaks and an absorption trough with width AA ~ 0.01 A in 
the rest frame. If this object is redshifted to z = 3, the width is widened to AA « 0.04 A in the 
optical band, which is beyond the observational limit of current instruments. 

(b) The surrounding medium subsequently forms a conical H II region, and a partially ionized 
envelope may be formed at the center. We assume that the scattering medium is hot (T = 10^ K) 
and the source is at the center of the slab. The emergent Lya profiles show symmetric double peaks 
and an U-shaped absorption trough at the line center. 

The trough has a width of AA 0.2 A, and the features that are formed in nearby starbursting 
galaxies are found to be observable using the current high resolution far UV spectroscopy such as 
HST/STIS. When those objects are located at z = 3, the width of the absorption feature is about 
1 A which can be marginally detected with the Keck LRIS. 

(c) In the next stage, an expanding superbubble develops. In this stage, there appear two Lya 
emission peaks. The major peak is formed by the central star cluster and the minor peak is formed 
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by the radiative and expanding supershell. Hence, the major peak is located at the hne center, and 
the minor peak is in the blue part of the line center. According to our investigation, the major and 
minor peaks may have an absorption trough at the center. 

The width of the trough corresponds to AA ~ 0.25 A. Hence, this central U-shaped absorption 
feature in nearby starbursting galaxies may be marginally observable using the HST/STIS. When 
those objects are located at z = 3, the width of the absorption feature becomes ~ 1 A, and this 
trough should be marginally detected by using ground based large telescopes. 

(d) When the neutralization of the supershell is complete, its H I column density reaches 
Nhi = 10^^~^^ cm~^, and its expansion gradually ceases. Hence, we modeled this stage by a hot 
source surrounded by a static and cold (T = 10 K) scattering medium. Since at this stage the 
central Lya source is extinguished, we considered the ultraviolet continuum source near Lya. The 
profiles emergent at this stage are also characterized by double peaks and U-shaped absorption 
trough. 

In this paper, we concentrated on the cases of moderate optical depths. However, the line 
center optical depth tq of the neutral envelopes in starbursting galaxies is expected to be 10^^^. In 
order for our model to be realistic, wc must investigate these extremely thick optical depths and 
include dust effects, which will enhance the width of absorption trough, and make its detection 
more feasible. We will study these cases in the future. 

In section 3.3 we show the angular distribution -P(/i) of escaping photons. In the case of a 
slightly thick medium, we get an almost isotropic distribution given by P(//) oc In a moderately 
thick medium, photons tend to escape more in the direction parallel to the slab and hence P{n) 
becomes convex upward. In the extremely thick medium the escaping photons arc emergent prefer- 
entially in the direction normal to the plane, and P{lJi) becomes convex downward. This means that 
when wing scatterings dominate the line transfer, the photons are affected by the geometrical shape 
of the scattering region and the quantum mechanical properties of scattering. Therefore we expect 
the development of polarization when Lya photons are radiatively transferred in an extremely thick 
aspherical medium (Chandrasekhar 1960; Lee k, Ahn 1998). This will be considered in the near 
future. 
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of Education. HWL is also supported by Brain Korea 21 awarded to Yonsei University. We thank 
Chris Pearson for careful reading of the manuscript. This work was supported by KOSEF Grant 
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Fig. 1. — Our emergent Lya profiles (solid lines) are compared with Neufeld's analytic solutions 
(dotted lines). The horizontal axis is the frequency in units of the thermal width, and the total 
flux of the line is normalized to I/Att in accordance with Neufeld's normalization. The profiles are 
symmetric about the origin, x = 0. Here the Voigt parameter is related with the temperature by 
a = 4.71 X 10~^T^ ^^'^ , where T^^ represents the temperature of medium in units of 10^ K. We note 
that the our results are in good agreement with analytical results by Neufeld for large qtq. 
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Fig. 2. — Last scattering positions of escaping photons. We can see the two peaks at Pz/tq 
and I-PzI/to ~ 0.95. See the text for the detailed explanation for the origin of these peaks. 
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Fig. 3. — Distributions of photon path lengths just before escape. The horizontal axis represents 

path lengths Tf, which is defined by the transverse optical depth in units of tq for the escaping 
photons in the slab after their last scattering. The vertical axis represents the probability density 
for a photon to have a last scattering path length in a range (t/,tj + dTf). There are three kinds 
of peaks: the first is at tj « 0, the second is at 1 < < 2, and the third is at tj « 2. 
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Fig. 4. — Angular distributions of emergent photons for various optical depths tq and Voigt pa- 
rameters a of the scattering medium. For isotropic distribution, P{n) oc /j,. We note that the 
curves become convex downward when otq is large, whereas they become convex upward when utq 
is smaller. 
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Fig. 5. — Distributions of the mean number of scatterings for emergent Lya photons. In the top 
panel (a), tq = 10^ and a = 4.3 x 10~^; in the bottom panel (b), tq = 10^ and a = 1.0 x 10~^. 




Fig. 6. — Emergent Lya profiles for the cases of a hot photon source surrounded by a cold medium, 
whose line center optical depths are shown in the figure. The Voigt parameter a = 1.49 x 10~^ 
corresponds to T = 10 K. Hence Ax = 1 corresponds to AA = 0.0016 A in the rest frame. The 
Lya profile of the source is assumed to be a Gaussian, exp(— x^/x^), where Xg = 30. We show the 
grobal profile in a small box in the figure. The total flux is normalized to be l/47r. The profiles 
are characterized by double peaks and the absorption trough around the center. The absorption 
trough ends at x ~ ±3, and hence the full width Ax ~ 6. The Lya profile in the small box is 
that for To = 10^, which shows a characteristic double peaks and the absorption trough at the line 
center. 
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Fig. 7. — Emergent Lya profiles for Case II, where monochromatic {xq = 0) sources he in the center 
and at and are surrounded by a hot medium. We assume that the slab-like scattering medium with 
the Voigt parameter a = 1.49 x 10~^ or the temperature T = 10 K typical of H II regions. We note 
that the profiles are symmetric to the origin. The central trough widens as the line center optical 
depth increases. Here Ax = 1 corresponds to AA = 0.05 A in the rest frame. 
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Fig. 8. — Emergent Lya profiles for sets of parameters of tlie scattering media. Here the monochro- 
matic Lya source is located at the center of the slab-like medium. In the figure we denote the Voigt 
parameters and line center optical depths of media by a and tq. The profiles are symmetric to 
X = and normalized so that the total flux becomes l/47r, and symmetric to x = 0. Note that 
when aro is large, wing scattering dominates the line transfer and results in broad bumps. Here 
Ax = 1 corresponds to AA = 0.05 A in the rest frame. 
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Fig. 9. — Emergent Lya profiles for Case III, in which monochromatic Lya sources are uniformly 
distributed in a hot slab-like medium. We assume here the Voigt parameter of the medium is 
a = 4.71 X 10^'^ which corresponds to T = 10^ K, which is a typical value for H II regions. The 
profiles arc symmetric to the origin. Comparing with the midplane case, we see that peaks become 
fatter to the line center. Note that here Ax = 1 corresponds to AA = 0.05 A in the rest frame. 
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Fig. 10. — Emergent profiles for various parameters of scattering media. Here the monochromatic 
Lya source is uniformly distributed in the slab-Hke medium. We denote in the figure the Voigt 
parameters and Hne center optical depths of media. The profiles are normalized so that the total 
flux is l/47r, and symmetric to a; = 0. Note that wing scatterings dominate the line transfer and 
result in broad bumps. The peaks are fatter to the line center that those of Case II, the midplane 
case. Here Ax = 1 corresponds to AA = 0.05 A in the rest frame. 
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Fig. 11. — Emergent Lya profiles for Case IV, in which the continuum source is located at the 
center of the slab-like hot medium with the Voigt parameter a = 4.3 x 10~^. We show the effects of 
the line center optical depths of the scattering media on the shapes of the Lya lines. The broader 
absorption troughs appear as tq grows larger. Here Ax = 1 corresponds to AA = 0.05 A in the rest 
frame, and therefore the width of the trough is about 0.3 — 0.5 A. 
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Fig. 12. — Emergent profiles for Case V, in which the continuum source is located at the center of 
the slab-like cold medium with temperature T = 10 K. We show the effects of the line center optical 
depths of the scattering media. As tq grows large, there appear the absorption troughs which are 
much broader than those in Case IV for the same tq. Here Ax = 1 corresponds to AA = 0.0016 A 
in the rest frame, and therefore the width of the trough ranges 0.01 - 0.03 A. 
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Table 1. Lya Emission Features of Starbursting Galaxies 



scatterer 



trough 



Stage 


Case 


To 


T 


AAo(A) 


AA3(A) 


profile type* 


(a) 


(I+V)3 


1Q3-5 


10 


0.01-0.03 


0.04-0.13 


T-I-DP 


(b) 


II+IV 


IQ3-5 


10" 


0.3-0.5 


1.2-2.0 


T-I-DP 


(c) 


(II+IV)i 










ST-hP Cyg 






103-5 


10-* 


0.25 


1 


T-hDP 


(d)5 


VI 


106-9 


101-2 






P Cyg 


(d) 


I?V3 


103-r. 


10 


0.01-0.03 


0.04-0.13 


T-I-DP 



Case I Hot source at the center surroun ded by cold scattering media 

Case II Monochromatic hot source at the center of a partially ionized hot medium 

Case III Distributed source in a slab w ith to < 10^ 

Case IV Continuum source at the center of a hot scatterer 

Case V Continuum source at the center o f a cold scatterer 

Case VI Scattered by expanding supershell (Future Work) 



1 Major emission line 

2 Minor emission line 

3 For reality tq = 10^'^^ should be achieved. 

* T = Absorption trough; DP = Double peaks; P Cyg = P Cyg type profile 
^ Future Work 



Table 2. Spectral Capabilities of Devices 



Instrument 


Spectrograph 


Range 


Resolution 






(A) 


(A) 


FUSEi 


HIST 


905-1187 


0.05-0.06 


HST/STIS2 


G140H 


1140-1740 


0.007-0.01 


Keck3 


LRIS 


3900-11000 


2-5 



iThe Fuse Observer's Guide Version 2.1 
2STIS Instrument Handbook v4.1 
^Oke et al. (1995) 



